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Abstract 

We apply two non-parametric methods to test further the hypothesis that log-periodicity 
characterizes the detrended price trajectory of large financial indices prior to financial crashes 
or strong corrections. The analysis using the so-called {H, (7)-derivative is applied to seven 
time series ending with the October 1987 crash, the October 1997 correction and the April 
2000 crash of the Dow Jones Industrial Average (DJIA), the Standard & Poor 500 and Nasdaq 
indices. The Hilbert transform is applied to two detrended price time series in terms of the 
\n{tc — t) variable, where tc is the time of the crash. Taking all results together, we find strong 
evidence for a universal fundamental log-frequency / = 1.02 ± 0.05 corresponding to the 
scaling ratio A = 2.67 ± 0.12. These values are in very good agreement with those obtained in 
past works with different parametric techniques. 

1 Introduction 



Several recent analyses [12, T^, 27, 28, 24, 23] have shown that large financial crashes are "outliers" 
or "kings" [22] with statistically different properties than the rest of the population of price returns. 
A proposed mechanism for understanding this "king" effect involves collective interactions between 
agents leading to a cascade of amplifications, which may result from several origins [^, 19, 20, 30, 
29], H, p. 



Associated with this "king" effect, log-periodic patterns [26] in financial price time series have 
been found to be precursory signatures of large crashes or large drawdowns [ ^ , ^, 33, 0, ^ 
m, |T5|, |, |17|, |] (see also ^ 



and references therein). Most of the analyses have been of a 
parametric nature, based on formulas involving combinations of powers laws and of log-periodic 
functions of the type cos[uj\n{tc 



t)] | |3T| , 13, 14], where tc is a "critical" time at or close to the 



crash. A first attempt at providing a non-parametric test of log-periodicity in financial time series 



was based on a spectral Lomb analysis [ ]25| ] of the oscillatory residuals to a preliminary power law 
fit to the time series [jl^, ^, 0. 

Here, we propose to extend the search for non-parametric signatures of log-periodicity. For 
this, we introduce two non-parametric techniques and apply them to seven financial time series 
culminating in a crash or a strong correction. The first non-parametric method is called the gen- 



eralized g-analysis or {H, g)-analysis and was introduced in Ref. [37] in the analogous context of 
the detection of log-periodic oscillations in precursory signals to material failures. The second 
non-parametric method uses the Hilbert transform to construct the local cumulative phase of the 
oscillatory component in the time series. By the new lights provided by these analyses that are cast 
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on the difficult empirical problem of detecting reliable precursory patterns to financial crashes or to 
strong corrections, we hope to provide what we believe are strong additional evidence for the reality 
of log-periodicity in price trajectories. Compared with a similar analysis of material failures using 



the {H, q)-analysis [p7[], we find that log-periodicity seems to be a significantly stronger signal in 



financial time series than it is for rupture. 

2 The generalized {H, -analysis 
2.1 Definitions 

The so-called g-analysis is based on the concept of a g-derivative, the inverse of the Jackson q- 



integral 011|], which is the natural tool ^ for describing discrete scale invariance []26|]. The 



g-derivative of an arbitrary function f{x) for any q G (0, 1) U (1, oo) is defined as 

where x is the distance to the critical point. If the signal is a price p{t) function to time t, then 
X = tc — t, where tc is the critical time, that is, the time of the culmination of the bubble, which is 



in general close to the time of the crash [ |20| , gjj. The case q —>■ I recovers the normal definition of 
derivative. It follows that 

Dyjix) = DJ{x/q). (2) 

It is thus enough to study Dqf{x) for q e (0, 1) to derive its values for all g's. 

The necessary and sufficient condition for a function f{x) to be homogeneous of order ip is 

Dqfix) = . (3) 

q — 1 x 

This expression suggests the introduction of a generalized g-derivative that we call {H, g) -derivative, 
such that the dependence in x of f{x) disappears for homogeneous functions, for the choice 
H = tp. Consider therefore the following definition 



such that D^~^f{x) recovers the standard g-derivative Dqf{x). For a power law function f{x) = 
Bx"^, DI^^"^[Bx"^] = B{1 — g™) /(I — g)™ is constant. For a statistically homogeneous function 

/(x) = Bx"^, D^^^f{x) = constant, where the symbol = means equivalence in distribution. 

The generalized {H, g) -derivative has two control parameters: the discrete scale factor q devised 
to characterize the log-periodic structure and the exponent H introduced to account for a possible 
power law dependence, i.e., to correct for the existence of trends in log-log plots. 

Let us consider the following log-periodic function of x = {tc — t)/T: 



A - Bt^x^ + Ct^x^ cos (u; In x) , (5) 



where the phase in the cosine = In r has been scaled away by a change of time scale. In the 
following, we shall refer to u) as the angular log-frequency and shall also use the log-frequency 
defined as 
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The {H, (7) -derivative of y{x) is 

D^y{x) = x^-"[B' + C'g{x)]. (7) 

where B' = -^f^^, C = and 

g{x) = cos{uj In x) — cos{ui In qx) (8) 
= Ci cos(w Inx) + C2 sin(L<j Inx), (9) 

where Ci = [1 — cos{uj In q)] > and C2 = sin(a;lnq'). In principle, we could optimize the 
value q to adapt it to the angular log-frequency lo to ensure Ci = 0, so that g{x) becomes a pure 
sinus: a study of the phase of the log-periodic oscillations as a function of q can in principle get 
access to lu. We have not tried this approach as it seems very sensitive to noise. 

Varying (3, we find that the analysis in terms of the [H, (7) -derivative of y{x) given by (js]) gives 
essentially the same results for the highest peaks {H, q) and their associated log-frequencies 
f{H, q) when varying (3 at fixed oj = 27: f but gives very different and volatile results when varying 
u; at fixed f3. This shows the value of the {H, (7)-analysis which is a robust detector of the log- 
frequency but not of the power law trend. 



2.2 Effect of errors in the critical time t, 



c 



An important issue is the determination of the critical time tc entering in the definition of the scaling 
parameter x = tc — t, such that the time series is (discretely) scale invariant with respect to scale 
transformation performed on x. Indeed, according to the theory of rational bubbles and crashes 



pOl , |2l| , |30| ], the crash is not a certain deterministic event but is described by its hazard rate h{t), 
which is the probability per unit time that the market crashes conditioned on the fact that it has not 
crashed yet. In this framework, the critical tc is the end of the bubble and coincides with the peak 
of the crash hazard rate. In practice, this means that the time of the crash, if it occurs, is expected 
to be close to and slightly smaller than tc- Thus, knowing the time of the crash does not provide 
a perfect measurement of tc- Therefore, there is a fundamental and unavoidable uncertainty in the 
determination of the value of tc entering the definition of x. This uncertainty is a severe limitation 



for parametric methods of detections of log-periodicity. In contrast, our previous work [ [37p has 
shown that the {H, (/) -derivative is rather insensitive to an error in the estimation of tc up to times 
of the order of the error made on the determination of tc (see figure 2 of ^%). Technically, this can 
be seen from the following derivation. 

An error in tc amounts to a spurious translation of the coordinate system {x, y) away from 
the origin (0,0) to a new origin (xo,yo)- This defined the following coordinate transformation 
(x,y) {x,y): 

y y ^ y = y - yo, 

where 7/0 = y{xo)- According to this transform, y{x) = y{x) — yo = y{x + xq) — yo and 
y{qx) = y{qx + xq) — yo = y{qx — qxo + xq). Thus, the {H, g) -derivative is 

j^H-/-\ Bq^ — B -\- C cos{u! In x) — Cq^ cos[u)ln{qx)] 
" ^^^^ " (1 - (7)^(x - xo)^x— ^ ^ 

where q = q + ii^i^^. 

Consider two values Xn+i > x„ of the control parameter x corresponding to successive local 
extrema of the log-periodic function, such that cos[aj ln(g„x„)] = cos[liJ ln(g„+ix„+i)], that is, 

ln{qnXn) - ln(^„+ix„+i) = 2k7r/uj . (12) 
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Therefore, 

Xn+l-qXo + XQ 2k2L ,,,, 

\ = e " (13) 

For Xn+i > Xn ^ Xq, wc obtain 

InXn+i —InXn = 2kiT/uJ , (14) 

that is, the translational error in the x variable has a negligible impact on the determination of the 
angular log-frequency u defined by ^ and recovered by the measure dlit), as long as the analysis 
is not too close to the critical point. In contrast, consider for instance xi < X2 comparable to xq 



and which satisfy Eq. (14) with = 1. In the translated coordinate system, the phase difference 
between these two points of the term Cq^ cos[lo ln{qx)] in the {H, g) -derivative of y{x) given by 



( |11D is given by 

A<^ = .;ln ^^-^^° + "V (15) 
xi - qxo + xo 

Since < g < 1, < A0 < 27r, which implies that the "period" near xq is stretched by the 
translation of x to the uncorrect critical value. In this case, the log-periodic oscillations disappear 
for X of the order of and less than xq (see figure 2 of [|^]). 

Suppose that the log-periodic oscillations for times tc — t < tc — tmax are not well-determined 
or unobservable. Then, for a given q, the {H, g)-derivative is not well-determined for 

t>t,- h^l^ . (16) 
q 

The smaller q is, the larger is the time interval in which the {H, g)-derivative is not determined. 
This explains our results below that small g's do not give good results and tend to over-emphasize 
slow log-periodic oscillations, hence the existence of the so-called and spurious B platform at very 
low log-frequency. The existence of spurious low log-frequencies is a general phenomenon which 
may derive from a variety of origins [^. 

In the sequel, we shall not attempt to optimize our analysis of tc and shall fix to the time 
of the crash or strong correction ending each bubble, keeping in mind that one should thus expect 
a deterioration of the searched signal close to the end of the time series. In practice, we have 
not found this to be a limitation for the present purpose of establishing the existence of the log- 
periodic patterns. However, the situation is not as bright when attempting to use this analysis in a 
prediction mode, that is, by scanning tc with an incomplete time series ending earlier than the crash. 



as discussed in [37]. 



2.3 Procedure and selection criteria 

Our analysis of a given time series consists in the following steps. 

1. Read the critical time tc as given by the time of the crash or of the strong correction and 
transform the time series p{t) into a function y{x) = p{t) with x = tc — t. 

2. Fix q and H to some arbitrary value and construct the corresponding {H, -derivative D^p{t) 
of the function y{x). 



3. Perform a spectral Lomb analysis [||, 25] of D^p{t) in order to identify the presence of 
possible periodic oscillations: the most significant oscillations occur at those log-frequencies 
that make the spectrum Pn{H, q) maximum. 

4. Change q and H and redo the analysis. In this work, we scan H E (—1, 1) with spacing 0.1 
and q E (0, 1) with spacing 0.05. 
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5. Having identified the strongest peak in the spectrum P/v {H, q) and its associated log-frequency, 
we study their dependence as a function of H and q. 



Let us illustrate our procedure on the daily evolution of the Dow Jones Industrial Average from 
02-Jan-1980 to the "Black Monday" on 19-Oct-1987, shown in Fig. [T[ 

In Figs. ^ and we present typical plots of the generalized g-derivative and its Lomb peri- 
odogram for the Dow Jones 1987 Crash. Figure ^ shows the evolution as a function of ln(tc — t) 
of the generalized g-derivative (left panel) and their Lomb periodograms (right panel) for fixed 
q = 0.65: (a-b) H = -0.9; (c-d) H = 0.1; and (e-f) H = 0.9. The log-periodic structures of 
Dqp{t) in plots (a), (c) and (e) are clearly visible. The amplitude of the log-periodic oscillations 
decreases with increasing H for fixed q. The envelops of the oscillations can be fitted well by the 
expression 

= I^B' ± C'^/c[+Cpj (17) 

giving the local maxima of the log-periodic function (js]). We observe a peak at /i = 1.06 in (b), 
fi = 1.06 in (d) and /i = 0.94 in (f). We also see harmonics /2 = 1.94 in (d) and /2 = 1.90 in 
(f). Figure | is for q = 0.9: (a-b) H = -0.9; (c-d) H = 0.1; and (e-f) H = 0.9. The Lomb peaks 
in plots (b) and (d) are very high with log-frequencies /2 = 2.02 and 1.96 that may be interpreted 
as the second harmonics 2/ of the frequency found in figure |[ We find a peak at /i = 1.15 in (d) 
and peaks at fi = 1.01, /2 = 1.87 and = 4.06, which makes more plausible the existence of 
harmonics of a fundamental log-periodic component in the signal. Figure ^ is for g = 0.3. The 
log-periodic structures of D^p{t) in plots (a), (c) and (e) are much less obvious visually than in 
previous figures. Nevertheless, the Lomb spectral analysis shown in panel (b) shows that the log- 
periodic structure in (a) is significant at two log-frequencies fi = 1.22 and /2 = 1.92. We note that 
the almost equal strength of the two overlapping spectral peaks at these two frequencies is probably 
causing a distortion in their estimation. 

In panels (d) and (f) of figure ^ the highest peaks occur at a significantly smaller log-frequency 
which is found close to 0.32. This log-frequency, which corresponds to approximately 1 to 1.5 
full log-periodic oscillation in the data, can usually be associated with noise rather than a genuine 
signal. Indeed, log-periodicity appears generically when dealing with power laws in the presence 



of some kind of cumulative noise [|9|| or with noise with long-range correlations p6|]. The most 
probable log-frequency resulting from noise is given by 

1.5 

ln(tc - tmin) - ln(tc - 

where tmin and tmax are respectively the times defining the beginning and the end of the time series. 
The value 0.32 of the highest peaks in panels (d) and (f) is compatible with /™p given by (H^). We 
should therefore exclude from our statistics those pairs {H, q) whose spectra are dominated by the 
spurious low log-frequency /™p resulting from noise. 

In addition, we should also discard from our statistics the pairs {H, q) whose log-frequencies / 
of the largest spectral peak are very large, because a large log-frequency corresponds to many os- 
cillations which most often can be associated with high-frequency noise. In Ref. [[T6|], log-periodic 
components with 27r/ > 14 were discarded for this reason. 

Finally, let us note that panels (d) and (f) of figure^ also contains secondary peaks at /i = 1.22 
and /2 = 1.92 in (d) and /i = 1.33 and /2 = 1.86 in (f) which are not far from the log-frequencies 
previously estimated. 
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3 Log-periodicity of Wall Street Crashes 



We now exploit the procedure described in section to analyze seven financial time series ending 
at 

1. the Dow Jones Industrial Average, October 1987 Crash, 

2. the Dow Jones Industrial Average, October 1997 strong correction, 

3. the S&P 500 Index, October 1987 crash, 

4. the S&P 500 Index, October 1997 strong correction, 

5. the Nasdaq Index, October 1987 crash, 

6. the Nasdaq Index, October 1997 strong correction, and 

7. the Nasdaq Index, April 2000 crash. 

In the spirit of R. Feynman |^, we now describe extensively the evidence we have accumulated 
on these seven time series, in particular by presenting all the evidence that we think relevant, so that 
the reader can judge for him- or her-self about the validity and robustness of the results. 

3.1 Time series ending at the Dow Jones October 1987 Crash 

The log-periodic patterns characterizing the bubble phase preceding this crash has been studied 



parametrically in [ |32| , |33| , |19| , |20| , gJJ. We synthesize the information obtained by our new non- 
parametric analyses such as those shown in Figs. ^, ^ and |^ for all pairs {H, q) by plotting in figure 
^ the log-frequency / of the largest spectral peak as a function of an index counting the pairs {H, q). 
Since we scan H G (—1, 1) with spacing 0.1 and q £ (0, 1) with spacing 0.05, this corresponds to 
324 pairs. 

• cluster Pi centered on /i = 1.04 ± 0.11 seems to be associated with a fundamental log- 
frequency in the data; 

• cluster P2 centered on /2 = 1.98±0.07 can be interpreted as the frequency 2/i of the second 
harmonics to the fundamental frequency /i ; 

• cluster B centered on fmp ~ 0.3 can be attributed to the spurious log-frequency defined in 



(18). 



Fig. ^presents a more detailed description of the log-frequency f{H, q) of the most significant 
peak in each Lomb periodogram of the {H, g)-derivative for the Dow Jones 1987 crash, by plotting 
/(i?, q) as a function of H and q. The clusters of figure § now corresponds to "platforms". It is 
convenient to classify all pairs of {H, q) into three classes based on the geometric shape of f{H, q): 
W (wedge or wall), P (platform) and B (bottom of valley or basin). In figure |^ we can clearly 
distinguish platforms Pi and P2 centered on /i = 1.04 it 0.11 and /2 = 1.98 ± 0.07. Two wedges 

^R. Feynman, "Surely you're joking, Mr Feynman", Norton, 1985, page 341: "If you're doing an experiment, you 
should report everything that you think might make it invalid - not only what you think is right about it: other causes that 
could possibly explain your results; ... Details that could throw doubt on your interpretation must be given, if you know 
them. You must do the best you can — if you know anything at all wrong, or possibly wrong — to explain it. If you 
make a theory, for example, and advertise it, or put it out, then you must also put down all the facts that disagree with it, 
as well as those that agree with it. ... In summary, the idea is try to give ALL of the information to help others to judge 
the value of your contribution; not just the information that leads to judgment in one particular direction or another." 
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are also indicated in Fig. |^ as Wi and W2 but their log-frequencies are not distinguishable from 
those of P2. The spurious cluster B is also clearly visible. 

In order to compare the relative values of different {H,q) pairs, It is also instructive to construct 
in Fig. ^ the spectral height PN{H,q) of the highest peak of each Lomb periodogram for each 
{H, (7) -derivative as a function of H and q. By matching Fig. ^ with Fig. 0, we verify that the two 
parts of B have quite high Lomb peaks, while Wi and W2 in between correspond to valleys in 
Fig. 0. In Fig. |^ pairs {H, q) with q > 0.5 lead to either fi or to its harmonic f2- Small values 
of q thus give a weaker signal for log-periodicity. This property will be seen to apply in all other 
time series discussed below. We also note that on average, for pairs {H, q) with q > 0.5, the largest 
Lomb spectral peak increases with decreasing H. This is also illustrated in Figs. ^ and 0. 

Define an optimal pair {H, q) satisfying that Pn{H, q) is the maximum of Pn{H, q). Figure 
^ shows on the left vertical coordinate and with the thin line with squares a sectional cut of Fig. ^ 
for H = —0.7, that is, plots the highest peak Pj\^{H,q) as a function of q. The right vertical 
coordinate and the thick line marked with circles correspond to a sectional cut of Fig. |^ for the 
same value H = —0.7 and gives the log-frequency of the highest peak as a function of q. The 
strongest spectral power is found for q = 0.6 for which the log-frequency is fi. 



3.2 Time series ending with the Dow Jones October 1997 strong correction 

The log-periodic patterns decorating the price trajectory of the bubble phase preceding this strong 
correction have been already studied parametrically in [^, 29]. 

Figure ^ shows the daily evolution of the Dow Jones Industrial Average from 02 -Jan- 19890 
to 19-Oct-1997. On Monday, September 27th, the Dow Jones index suffered from a strong loss 
of more than 7% that ended a strong bull regime of strongly appreciating prices. Many observers 
believed on this day that it was the starting point of a crash. The next day, a rally of +5% alleviated 
fears and started a market phase that last approximately 3 months in which the market remained 
essentially flat. The results of the analysis of the {H, g)-derivative of this time series shown in 



figures 10 to 16 are found to be very similar to those reported previously for the Dow Jones time 



series of ending with the October 1987 crash. Two novel features appear. First, figure |13| shows the 
existence of an additional cluster P1/2 at fi/2 = 0.59 ± 0.04, that it is tempting to interpret as the 
first sub-harmonic of the fundamental log-frequency /i. Second, for q = 0.9, we observe in figure 



1 1 several harmonics beyond the second harmonics, suggesting a very strong log-periodicity with 



more complex nonlinearity that for the October 1987 crash. 



3.3 Time series ending with the S&P 500 Index, October 1987 crash 



Figure |17| shows the Daily evolution of the S&P 500 Index from 02-Jan-I980 to the "Black Mon- 
day" on 19-Oct-1987, which started the crash of October 1987. The time series is obviously not 
independent of the Dow Jones time series shown in figure |l| as the S&P 500 Index contains the 30 
blue chip stocks constituting the Dow Jones index. However, it is influenced by the behavior of 
470 other stocks and we think it useful to carry out a detailed analysis of this time series as well, to 
compare with the Dow Jones time series and develop a sense about how different realization of the 
noise and of the signal may modify the log-periodicity. 



The results of the analysis of the {H, g) -derivative of this time series shown in figures 18 to 23 
are found to be very similar to those reported previously in the two previous time series. The same 
fundamental log-frequency /i and its harmonics are observed. The only really novel feature is the 
appearance of what is probably a spurious "wedge" cluster W with log-frequency fw = 1.74ib0.07 



shown in figure 20 
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3.4 Time series ending with the S&P 500 Index, October 1997 strong correction 

Figure || shows the daily evolution of the S&P 500 Index from 02-Jan-1990 to 27-Oct-1997. Sim- 
ilarly to the comparison of the Dow Jones time series ending in October 1987 and 1997, we think it 
useful to compare the S&P 500 time series ending in October 1987 and 1997, to test for the robust- 
ness of the log-periodic signal with respect to modifications and/or alterations brought by different 



realization of the noise and of the signal. Figures ^ to g9| gives the results of the analysis of the 
{H, (/) -derivative of this time series. The results are highly consistent with those obtained in the 
previous three time series. 

3.5 Time series ending with the Nasdaq Index, October 1987 crash 



Figure 30 shows the daily evolution of the Nasdaq Index from ll-Oct-1984 to the "Black Monday" 



on 19-Oct-1987, which started the crash. Figures 31 to 36 show the results of our analysis, which 



essentially recovers similar results as for the previous time series. The main difference is the weaker 
evidence of clustering into three clusters Pi with /i = 0.88 ± 0.04, P2 with /2 = 1.45 ± 0.06 and 
B with fs ~ 0.21 (which corresponds as in other previously analyzed time series to approximately 
one oscillation). 

3.6 Time series ending with the Nasdaq Index, October 1997 strong correction 



Figure 37 shows the daily evolution of the Nasdaq Index from 02-Jan-1990 to 27-Oct-1997, at 



which a strong correction occuned. Figures |38| to ^ show the results of our analysis of this time 



series, which are very similar to previous ones. 

3.7 Time series ending with the Nasdaq Index, April 2000 crash 

Figure^ shows the time evolution of the daily quotation of the Nasdaq Index from 03-Mar-1997 
to 14-Apr-2000. This crash corresponds to the burst of the so-called "new economy" bubbles that 



has spurred in the few preceding years. We refer to [|15|, g9p for a detailed discussion of this crash 
and a thorough analysis of its log-periodic signature using a parametric approach as well as the 
non-parametric Lomb spectral analysis of the residuals to the fit of the Nasdaq index time series to 
a power law. Figures to ^ show the results of our analysis, which are very similar to previous 
time series. 



4 Summary and discussion of the {H, -analysis 

By comparing the log-frequency f{H, q) associated with the highest peak Pn{H, q) of the Lomb 
spectrum of the {H, g) -derivative for the seven financial indices reported above, we find striking 
similarities, in particular in the shape and values of the plateaux and clusters of log-frequencies. 
Some of these similarities can be attributed to the fact that all these time series are not indepen- 
dent, especially those covering the same time spans which use indices with overlapping stocks 0. 
However, the found universality of the fundamental log-frequency is quite remarkable as it seems 

■^The DJIA is an index of 30 "blue-chip" US stocks. It is tlie oldest continuing US market index. It is called an 
"average" because it was originally computed by adding up stock prices and dividing by the number of stocks (the very 
first average price of industrial stocks, on May 26, 1896, was 40.94) and should ideally represent a correct measure of the 
state of the economy. The methodology remains the same today, but the divisor has been changed to preserve historical 
continuity. The standard & Poor 500 index is a capitalization weighted average of 500 of the major stocks on the US 
market. Thus the S&P 500 includes the 30 stocks constituting the Dow Jones index and it is natural to expect strong 
correlations between the two, especially at times of strong herding as in speculative bubbles. 
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Table 1 : Synthesis of the results of the analysis of seven price time series using the {H, g) -derivative. 
The periods of the analyzed data sets are: from 02-Jan-1980 to 19-Oct-1987 for the Dow Jones and 
S&P 500 indices, from ll-Oct-1984 to 19-Oct-1987 for the Nasdaq index, from 02-Jan-1990 to 
27-Oct-1997 for the Dow Jones, S&P 500 and Nasdaq Indices, and from 03-Mar-1997 to 14-Apr- 
2000 for the Nasdaq index. The columns /i, /2, /s and fi/2 are respectively the fundamental 
log-frequency, its second and third harmonics and its first sub-harmonic. The columns {H,q), f 
and Pn are the optimal pairs and their corresponding log-frequency and their Lomb peak height. 



independent of the specific realization of the noise decorating the log-periodic patterns, as proxies 
by the different indices for the same time periods, and is robust over the three time periods that have 
been investigated. 

We summarize in Table |^ all the results presented in previous figures by listing the most signif- 
icant log-frequencies for each of the seven time series. Apart from the Nasdaq time series ending 
on Oct. 1987, all the other 6 time series give a "fundamental" log-frequency which is remarkably 
universal and compatible with a universal value close to 1. This value corresponds to a preferred 
scaling ratio between successive maxima of the oscillations equal to A = e^/-^ 2.7. One can also 
observe the presence of harmonics and of the first sub-harmonic of this fundamental log-frequency, 
which reinforce the evidence for the existence of a discrete scaling invariance: indeed, the sub- 
harmonic corresponds to A^, that is, to two steps of the discrete hierarchy; as for the harmonics, any 
nonlineaiity in the price dynamics will produce them from a fundamental log-frequency. Taking 
into account of the log-frequencies /i, /2 and /a, excluding the Nasdaq 1987 time series, and as- 
suming that all these log-frequencies derive from a universal fundamental value /, our best estimate 
for / is then / = 1.02 ± 0.05. The corresponding preferred scaling ratio is then A = 2.67 ± 0.12, 



which is in excellent agreement with that from other analyses [ p3i [27p . The standard deviation of 
A is evaluated by the formulae a\ = e^l^ af/f"^. Our analysis also suggest that the most sensitive 
range for q lies in the interval q G [0.6, 0.85]. 

The Nasdaq 87 time series seems to belong to a different class with a significantly smaller 
fundamental frequency /i. This may be due to the fact that the Nasdaq was more like an "emerging" 
market before the October 1987 Crash, and thus shared the property of having larger fluctuations, 
as found in a recent systematic study of other emerging markets [p^]. 

5 Analysis of log-periodicity using the Hilbert transform 
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5.1 Definition and useful properties of the Hilbert transform 

The Hilbert transform (HT) of a signal x{t) is defined by the equation 



Hs) = - , (19) 

vr J-oo s — t 



where the integral is the Cauchy principal value. The reconstruction formula 



a:(t) = -- ^ds (20) 
vr t — s 



defines the inverse HT. Expression (19) shows that x{s) = (1/vrt) ® x{t) is the convolution of 



(l/vrt) with x{t), which by taking the Fourier transform implies that 

X{u) = -j sign{u)X{iy) , (21) 

where = —1 and X{ij) (respectively X(zv)) is the Fourier transform of x{t) (respectively x(s)). 
This means that the HT produces a — 7r/2 radian phase shift for the positive frequency components 
of the input x{t), while the amplitude does not change. A signal x{t) and its HT x{s) have the same 
amplitude spectrum, the same autocorrelation function, they are orthogonal and the HT of x(s) is 
-x{t). 

The HT is thus an ideal tool for detecting an arbitrary phase function P{t) in a signal of the 
form x{t) = A{t) cos{6{t)) + n{t), where n{t) is a noise with small amplitude (n <C A), and A{t) 
is a slowly varying amplitude modulation. The phase 9{t) is in principle recovered by the formula 

9{t) =tan-^[x{t)/x{t)] . (22) 

However, the phase obtained in this way consists in saw-tooth peaks that drop sharply from tt down 
to — vr. It is thus necessary to unwrap the phase to obtain a smooth continuous phase. 



An example is given in figure |49| which shows the function 



2 sin(27rx) - 4, < x < 3, 
x{t) = e(t) + { sin(47rx), 3 < x < 10, (23) 
3sin(27rx), 10 < x < 20. 

where e(t) is normally distributed, which consists of a high-frequency burst in a low-frequency 



background. The unwrapped phase 9{t) is plotted in Fig. 50. The abrupt translation of x{t) at 
t = 0.3 leads to an incorrect determination of the phase in this interval, and the corresponding 
distortion also spreads in the high-frequency interval approximately up to t = 5.5. Otherwise, the 
abrupt change of frequency is correctly detected at i = 10 notwithstanding the abrupt change in 
amplitude. 

5.2 Hilbert transform to financial time series 

Before applying the HT transform, we prepare the data (financial time series) by transforming the 
time variable t into x = ln(tc — t), where tc is fixed as in the {H, g)-analysis to the time of the 
realized crash. The new variable r is such that log-periodic oscillations in the variable t become 
regular periodic oscillations in the variable x. We then fit the price time series as a linear function of 
X. The residuals of this fit for the two S&P 500 price time series ending with the October 1987 and 
the October 1997 correction are shown in figures ^ and We then apply the Hilbert transform 



to these residuals and obtain the corresponding unwrapped phase in figures 53 and 54. Apart from 



jumps and distortions at the boundaries due to fast oscillations and rapid changes of amplitudes, the 
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unwrapped phases provide a reliable determination of the log-frequency in the vicinity of / = 1, in 
agreement with the results obtained above with the {H, (7)-analysis. 

As an attempt to improve the HT, we first smooth the financial time series in two steps, by first 
performing a spline interpolation with 100 evenly spaced sampling points in ln(tc — t) and then 
applying the Savitzky-Golay filter with centered window of width 21 and six-order polynomials. 
Such smoothing will alleviate some of the problems in the HT coming from fast-varying amplitudes. 
The resulting residuals (that is, filtered versions of the curves shown in figures |1| and 52 are given 



in figures |55| and q6[ The corresponding unwrapped phases obtained by taking the HT of these data 
are shown in figures ^ and The slope of these unwrapped phases provide a good estimation of 
the log-periodic frequency close to 1, in very good agreement with the {H, q) analysis. 



6 Conclusion 

In summary, we have applied two non-parametric analyses to test the hypothesis that log-periodicity 
is present in financial time series preceding crash or strong corrections. The {H, -derivative and 
the Hilbert transform performed on the detrended price in the ln{tc — t) variable are both consistent 
with the existence of a strong log-periodic signal with log-frequency / ~ 1. 

Acknowledgments: We are grateful to Vladilen Pisarenko for constructive help in Hilbert 
transform. This work was partially supported by the James S. Mc Donnell Foundation 21st century 
scientist award/studying complex system. 
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Figure 1: The daily evolution of the Dow Jones Industrial Average from 02- Jan- 1980 to the "Black 
Monday" on 19-Oct-1987. 
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Figure 2: Time dependence in the variable ln(tc — t) of the {H, g)-derivative (left panels) for the 
Dow Jones time series shown in figure [l] and their Lomb spectra (right panels) for fixed q = 0.65 
and varying H: (a-b) H = —0.9; (c-d) H = 0.1; and (e-f) H = 0.9. The log-periodic structures 
of D^p{t) in plots (a), (c) and (e) is visible to the naked eye. The amplitude of the log-periodic 
oscillations decreases with increasing H for fixed q. We observe peaks at /i = 1.06 in (b), at 
/i = 1.06 in (d) and at /i = 0.94 in (f). We also find harmonics /2 = 1.94 in (d) and /2 = 1.90 in 
(f). Recall the definition (||,^ of the log-frequency. 
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Figure 3: Same as Fig. | but for q = 0.9: (a-b) H = -0.9; (c-d) H = 0.1; and (e-f) H = 0.9. The 
peaks of the Lomb spectra in plots (b) and (d) are very high with log-frequencies /2 = 2.02 and 
1.96 that can be interpreted as harmonics 2/ of the log-frequencies observed in Fig. ^. We also see 
a peak at /i = 1.15 in (d) and peaks at /i = 1.01, /2 = 1.87 and = 4.06 in (f), suggesting the 
observation of several levels a discrete hierarchy. 
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Figure 4: Same as Fig. ^ but for q = 0.3. The log-periodic structures of D^p{t) in plots (a), (c) and 
(e) are much weaker than in previous figures. Plot (b) illustrates that the log-periodic structure in 
(a) is nevertheless very significant with the log-frequencies /i = 1.22 and /2 = 1.92 corresponding 
to the highest peaks of the spectrum. In the plots (d) and (f), we also observe secondary peaks at 
/i = 1.22 and /2 = 1.92 in (d) and /i = 1.33 and /2 = 1.86 in (f) that are consistent with previous 
values. However, the highest peak is in both cases found for a very low log-frequency /'"^ = 0.32, 
that can be shown to result from the trend induced by this choice of q. See text for an explanation. 
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Figure 5: For the Dow Jones time series shown in figure we report all the log-frequencies corre- 
sponding to the highest peak of the Lomb spectrum obtained for each of the pair {H, q) that have 
been systematically sampled with g] — 1, 1[ with spacing 0.1 and with q g]0, 1[ with spac- 
ing 0.05, providing a total of 18 x 18 "best" log-frequencies. The abscissa is counting the 324 
pairs with an arbitrary choice of indexing. Three clusters are clearly delineated: cluster Pi with 
/i = 1.04 it 0.11, cluster P2 with /2 = 1.98 ± 0.07 and cluster B with a log-frequency close to 
0.3. 
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Figure 6: Dependence of the log-frequency f{H,q) of the most significant peak in each Lomb 
spectrum of the {H, g) -derivative for the Dow Jones financial time series ending at the October 
1987 crash. Two platforms Pi and P2 with /i = 1.04 ± 0.11 and /2 = 1.98 ± 0.07 can be 
observed. 
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Figure 7: Dependence of the highest peak Pj\[{H, q) in each Lomb spectrum of the {H, g) -derivative 
for the Dow Jones time series ending at the October 1987 crash. 
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Figure 8: Right vertical axis and thin Une marked with squares: dependence as a function of q for 

fixed H = —0.7 of the highest peak Pn{H, q) of the {H, g') -derivative of the Dow Jones time series 
ending at the October 1987 crash. Left vertical axis and thick line marked with circles: dependence 
as a function of q for fixed H = —0.7 of the log-frequency associated with these highest spectral 
peaks. 
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Figure 9: The daily evolution of the Dow Jones Industrial Average from 02-Jan-1990 to 27-Oct- 
1997. 



22 




Figure 10: The evolution of the generalized g-derivative (left panels) for the Dow Jones time series 
ending at the October 1997 strong correction and their Lomb spectra (right panels) for fixed q = 0.7: 
(a-b) H = -0.9; (c-d) H = 0.1; and (e-f) H = 0.9. The highest peaks in the right panels are: (b) 
/i = 0.96 and /a = 2.34; (d) /i = 1.05 and /s = 2.31; and (f) /i = 1.04 and /s = 2.21. 
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Figure 11: Same as Fig. [ig for g = 0.9: (a-b) H = -0.9; (c-d) H = 0.1; and (e-f) H = 0.9. The 
log-frequencies and the associated harmonics are remarkably well-defined in the Lomb spectra: (b) 
/i = 1.10, /2 = 2.36, /s = 3.18, h = 4.21, h = 5.74, /g = 6.64 and fj = 7.21; (d) /i = 1.12, 
/2 = 2.35, /3 = 3.21, U = 4.28, h = 5.19, h = 5.87 and /r = 7.22; and (f) /i = 1.10, 
/2 = 2.25, /a = 3.19 and fi = 3.91. The amplitude of the log-periodic oscillations decreases with 
increasing H for fixed q. 
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Figure 12: Same as Fig. [Tg for g = 0.4: (a-b) H = -0.9; (c-d) H = 0.1; and (e-f) H = 0.9. The 
log-periodic oscillations are ambiguous and the extracted log-frequencies in (d) and (f) are close to 
the spurious log-frequency f^^. We can also see peaks at fi/2 = 0.52 in (b), /^/2 = 0-61 in (d) 
and fi/2 = 0.65 in (f), which may be probably interpreted as sub-harmonics of the fundamental 
log-frequency /i. 
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Figure 13: Same as figure ^ for tlie Dow Jones time series ending at the October 1997 strong 
correction. Four clusters Pi with /i = 1.04 it 0.05, P2 with /2 = 2.34, P1/2 with ^^12 = 
0.59 lb 0.04 and B are clearly visible. The small log-frequency = 0.21 of cluster B corresponds 
to only one oscillation in the g) -derivative. 
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Figure 14: Dependence of the log-frequency f{H,q) of the most significant peak in each Lomb 
periodogram of the [H, (7) -derivative for the Dow Jones time series ending with the October 1997 
strong correction. The two platforms Pi and P2 have log-frequencies of /i = 1.04 it 0.05 and its 
2/ harmonic /2 = 2.34, while P1/2 shows probably what can be interpreted as the subharmonic 
/1/2 = 0.59 ± 0.04. 
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Figure 15: Dependence of the highest peak PN{H,q) of each Lomb periodogram of the {H,q)- 
derivative for the Dow Jones time series ending with the October 1997 strong correction. 
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Figure 16: Dependence as a function of q for fixed H = —0.8 of the highest spectral peak Pn{H, q) 
of the (iJ, g) -derivative, shown on the left vertical axis with thin line marked with squares and of 
the associated log-frequencies f{H, q) on the right vertical axis with the thick Une marked with 
circles, for the Dow Jones time series ending with the October 1997 strong correction. 
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Figure 17: Daily evolution of the S&P 500 Index from 02-Jan-1980 to the "Black Monday" on 
19-Oct-1987. 
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Figure 18: Left panels: evolution of the generalized q-derivative for the S&P 500 time series ending 
with the October 1987 Crash and their Lomb periodograms (right panel) for fixed q = 0.6: (a-b) 
H = -0.9; (c-d) H = 0.1; and (e-f) H = 0.9. The log-periodic structures of D^p{t) in plots 
(a), (c) and (e) are clearly visible. The amplitude of the log-periodic oscillations decreases with 



increasing H for fixed q. The envelops of the oscillations can be fitted very well with Eq. (17). 
We see Lomb peaks at /i = 1.03 in (b), /i = 1.06 in (d) and /i = 0.91 in (f) which are very 
significant. We can also observe the harmonics /2 = 1.65 in (b) which is reddened by the highest 
peak and /2 = 1.92 in (d). 
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Figure 19: Same as Fig. for g = 0.85: (a-b) H = -0.9; (c-d) H = 0.1; and (e-f) H = 0.9. 
The spectral peaks in panels (b) and (d) are very high with log-frequencies /2 = 2.03 and 1.97 
corresponding to the second harmonics of /i. We can also see a peak at /i = 1.14 in (d) and peaks 
at /i = 0.99, /2 = 1.90, /s = 3.10 and fi = 4.10 in (f) that suggest a strong log-periodicity with 
its several harmonics. 
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Figure 20: Same as figure ^ for tiie S&P 500 time series ending at the October i987 crasli. Four 
clusters can be observed: Pi witli /i = 1.00 ± 0.11, P2 with /2 = 1.97 ± 0.06 and B with 
Jb = 0.21 probably associated with the spurious log-frequency fmp', the fourth "wedge" cluster 
W has log-frequency fw = 1-74 ± 0.07. 
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Figure 21: Dependence of the log-frequency f{H, q) of the most significant peak in each Lomb 
spectrum of each {H, (7)-derivative for the S&P 500 time series ending with the October 1987 crash. 
We observe two platforms Pi and P2 with /i = 1.00 ± 0.11 and /2 = 1.97 ± 0.06 respectively. 
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Figure 22: Dependence of the highest peak Pj^{H,q) in each Lomb spectrum of the {H,q)- 
derivative for the S&P 500 time series ending with the October 1987 crash. 
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Figure 23: Dependence as a function of q for fixed H = —0.3 of the highest peak Pn{H, q) shown 
as the thin line marked with squares with scales on the left vertical axis and of the associated log- 
frequencies f{H, q) shown as the thick line marked with circles with scales on the right vertical 
axis, for the S&P 500 time series ending with the October 1987 crash. 
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Figure 24: Daily evolution of the S&P 500 Index from 02- Jan- 1990 to 27-Oct-1997. 
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Figure 25: Evolution of the {H, (7) -derivative (left panel) for the S&P 500 time series ending with 
the 1997 strong correction and their Lomb spectra (right panel) for fixed q = 0.65: (a-b) H = —0.9; 
(c-d) H = 0.1; and (e-f) H = 0.9. The log-periodic structures of D^p{t) in plots (a), (c) and (e) 
are clearly visible. The amplitude of the log-periodic oscillations decreases with increasing H for 
fixed q. The envelops of the oscillations can be fitted very well with Eq. (|l7]). We see Lomb peaks 
at /i = 1.07 in (b), /i = 1.06 in (d) and /i = 1.06 in (f), which are very significant. We can also 
distinguish the sub-harmonics in the right panels. 
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Figure 26: Same as figure ^ for tlie S&P 500 time series ending at tiie October 1997 strong correc- 
tion. Three clusters can be observed: Pi with /i = 1.06 ± 0.04, P1/2 with /1/2 = 0.61 ± 0.05 
and B with Jb = 0.21. 




Figure 27: Dependence of the log-frequency f{H, q) of the most significant peak in each Lomb 
periodogram of the {H, -derivative for the S&P 500 time series ending with the October 1997 
strong correction. We see clearly two platforms Pi and P1/2 with /i = 1.06 ± 0.04 and = 
0.61 ±0.05. 
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Figure 28: Dependence of highest peak Pj^{H,q) in each Lomb periodogram of the {H,q)- 
derivative for the S&P 500 time series ending with the October 1997 strong correction. 
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Figure 29: Dependence with qiov H = —0.5 of the highest spectral peak Pn{H, q) shown as the 
thin line marked with squares with left vertical scales and of the associated log-frequencies f{H, q) 
shown as the thick line marked with circles with right vertical axis, for the S&P 500 time series 
ending with the October 1997 strong correction. 
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Figure 30: Daily evolution of the Nasdaq Index from ll-Oct-1984 to the "Black Monday" on 19- 
Oct-1987. 
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Figure 3 1 : The time evolution of the generalized g-derivative (left panels) for the Nasdaq index time 
series ending with the October 1987 Crash and their Lomb spectra (right panels) for fixed q = 0.5: 
(a-b) H = -0.9; (c-d) H = 0.1; and (e-f) H = 0.9. Three log-periodic oscillations of D^p{t) 
in plots (c) and (e) are clearly visible, while (a) exhibits a strong trend. The amplitude of the log- 
periodic oscillations decreases with increasing H for fixed q. The envelops of the oscillations can be 



fitted very well with Eq. (|17[). In the right panels, we observe peaks at the following log-frequencies: 
/i = 0.85 in (b), /i = 0.88 in (d) and /i = 0.84 in (f). 
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Figure 32: Same as Fig. |3l| for g = 0.9: (a-b) H = -0.9; (c-d) H = 0.1; and (e-f) H = 0.9. The 
log-periodic structures of D^p{t) in plots (a), (c) and (e) are clearly visible. The two highest Lomb 
peaks lie at /a = 1.60 in (b), /a = 1.39 in (d) and /a = 1.37 in (f). 
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Figure 33: Same as figure || for tiie Nasdaq time series ending at the October 1987 crasli. Three 
clusters can be observed: Pi with /i = 0.88 ± 0.04, P2 with /2 = 1.45 ± 0.06 and B with 
Jb ~ 0.21 (which corresponds as in other previously analyzed time series to approximately one 
oscillation). Compared to previous time series, the clusters are less clearly defined with a continuous 
range of log-frequencies joining clusters Pi and P2. 
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Figure 34: Dependence of the log-frequency f{H, q) of the most significant peak in each Lomb 
periodogram of the {H, g)-derivative for the Nasdaq index time series ending at the October 1987 
crash. One can observe two types of platforms Pi and P2 with /i = 0.88 ± 0.04 and /2 = 
1.45 ±0.06. 
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Figure 35: Dependence of the highest peak Pj^{H, q) in each Lomb periodogram of the {H, q)- 
derivative for the Nasdaq index time series ending at the October 1987 crash. 
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Figure 36: Dependence wit q for fixed H = 0.4 of the highest peak Pn{H, q) of the Lomb spectrum 
(thin line marked with squares) with scales given by the left vertical axis and of the associated log- 
frequencies f{H, q) shown as the thick line marked with circles with scales given by the right 
vertical axis, for the Nasdaq index time series ending with the October 1987 crash. 
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Figure 37: Daily evolution of the Nasdaq Index from 02-Jan-1990 to 27-Oct-1997. 
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Figure 38: Time evolution of the generalized g-derivative (left panels) for the Nasdaq index ending 
with the October 1997 correction and their Lomb periodograms (right panels) for fixed q = 0.65: 
(a-b) H = -0.9; (c-d) H = 0.1; and (e-f) H = 0.9. The Lomb peaks at fi = 1.14 in (b), /i = 1.13 
in (d) and /i = 1.10 in (f) are very significant. 
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Figure 39: Same as figure ^ for tfie Nasdaq time series ending at the October 1997 correction. 
Several clusters can be observed: Pi with log-frequency fi = 1.11 ± 0.04, P1/2 with /1/2 = 
0.51 ± 0.07, W P3 with /g = 2.90 ± 0.03 and B with Jb = 0.21. See figure || for a better 
definition of the clusters. 
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Figure 40: Dependence of the log-frequency f{H, q) of the most significant peak in each Lomb 
periodogram of the {H, g) -derivative for the Nasdaq index time series ending with the October 
1997 crash. One can identify several structures: platform Pi with log-frequency /i = 1.11 it 0.04, 
Pi/2 with fi/2 = 0.51 ± 0.07, W + Ps with /g = 2.90 ± 0.03. 
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Figure 42: Dependence with q for fixed H = —0.4 of tiie highest peak Pj\[{H, q) shown as a thin 
line marked with squares with left vertical axis and of the associated log-frequencies f{H, q) shown 
as the thick Une marked with circles with right vertical axis, for the Nasdaq index time series ending 
at the October 1997 correction. 
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Figure 43: Daily evolution of the Nasdaq Index from 03-Mar-1997 to 14-Apr-2000. 
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Figure 44: Time evolution of the generalized g-derivative (left panels) for the Nasdaq index time 
series ending with the April 2000 crash and their Lomb periodograms (right panels) for fixed q = 
0.7: (a-b) H = -0.9; (c-d) H = 0.1; and (e-f) H = 0.9. 
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Figure 45: Same as figure ^for the Nasdaq time series ending at the April 2000 crash. Several 
clusters can be observed: Pi with log-frequencies /i = 0.98 ±0.04, P1/2 with/1/2 = 0.48 ±0.13 
and B with fs = 0.25. Two small clusters are also observed whose significance is not obvious. 
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Figure 46: Dependence of the log-frequency f{H, q) of the most significant peak in each Lomb 
periodogram of the {H, (7)-derivative for the Nasdaq index time series ending with the April 2000 
crash. The two main plateaux are Pi with log-frequency /i = 0.98 ± 0.04, P1/2 with /1/2 = 
0.48 ±0.13. 
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Figure 47: Dependence of the highest peak Pj^{H, q) in each Lomb periodogram of the {H,q)- 
derivative for the Nasdaq index time series ending with the April 2000 crash. 
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Figure 48: Dependence with q for fixed H = —0.5 of tlie liigliest peak Pn{H, q) sliown as tlie thin 
line marked with squares with the left vertical axis and of the associated log-frequencies f{H, q) 
shown as the thick line marked with circles with the right vertical axis, for the Nasdaq index time 
series ending with the April 2000 crash. 
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Figure 51: Difference between the S&P 500 price time series ending with the October 1987 crash 
and the fit with A — B ln(tc — t), where A and B are two constants, as a function of ln(tc — t). tc 
is the time of the crash. 
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Figure 52: Same as figure gjj for the S&P 500 price time series ending with the October 1997 
correction. 
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Figure 53: Unwrapped phase 6 extracted from the residuals in Fig. |1] using the Hilbert transform. 
The slopes of the three linear segments are respectively the log-frequencies / = 0.9, / = 0.9 
and / = 1.9. This is in agreement with the existence of a fundamental log-frequency and its 
harmonics. The jumps in the phase are due to the fast oscillations of small amplitudes in the interval 

Htc-t) G (5,7). 




Figure 54: Unwrapped phase 6 extracted from the residuals shown in Fig. 52 using the Hilbert 
transform. Two parallel linear segments are indicated by the straight lines whose slopes « 1.0 
give the value of the log-frequency. The jumps are caused by the fast fluctuations in the second 
oscillation in Fig. 52 near ln(tc — t) = 5.7. The fluctuations in the fourth oscillation for ln(tc — > 
7.4 in Fig. completely spoil the Unear dependence of the unwrapped phase. 



60 




;l , , , , , , , 1 

4 4.5 5 5.5 6 6.5 7 7.5 8 



Figure 55: The residuals for the October 1987 crash time of the S&P 500 index from the smoothened 
data using spline interpolation and then Savitzky-Golay filter. The 100 interpolation points are 
evenly sampled in ln{tc — t). 
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Figure 56: Same as Fig. ^ but for the October 1997 correction of the S&P 500 index. 
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Figure 57: The unwrapped phase 6 extracted from the residuals in Fig. q5| using the Hilbert trans- 
form. Small oscillations within one oscillation lead to the undulations in the unwrapped phase. The 
two straight lines define log-frequencies equal respectively to / = 0.85 and / = 1.8. 
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Figure 58: The unwrapped phase extracted from the residuals in Fig. 56 using the Hilbert trans- 
form. The main part of the phase can be fitted to a straight line of slope giving the log-periodic 
frequency / = 1.0. Small oscillations within one oscillation lead to the undulations in the un- 
wrapped phase. 
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